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We theoretically study the effect of atomic nonlinearity on the tunneling time in the case of an 
atomic Bose-Einstein condensate (BEC) traversing the laser-induced potential barrier. The atomic 
nonlinearity is controlled to appear only in the region of the barrier by employing the Feshbach 
resonance technique to tune interatomic interaction in the tunneling process. Numerical simulation 
shows that the atomic nonlinear effect dramatically changes the tunneling behavior of the BEC 
matter wave packet, and results in the violation of Hartman effect and the occurrence of negative 
tunneling time. 
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Quantum tunneling of a wave packet through a poten- 
tial barrier is one of the fundamental topic in quantum 
physics (Tj-Q . The issue of tunneling time has attracted a 
lot of attention for decades since it was first put forward 
by Condon [lj . In the early 1960's, Hartman predicted 
that the tunneling time becomes independent of barrier 
length for thick enough barriers, ultimately resulting in 
unbounded tunneling velocities [5| ■ Such a phenomenon, 
termed as the Hartman effect later on, seems to imply 
the superluminal velocities inside the barriers and leads 
to a wide interest in many different fields 

Mathematically, the quantum tunneling is governed by 
the Schrodinger equation, which is a linear equation de- 
scribing the quantum wave nature of a single particle. So 
far the Hartman effect or related topics studied in the lit- 
erature are limited to the single-particle linear case. Af- 
ter mid-1990's, there are significant advancements in the 
realization of atomic Bose-Einstein condensate (BEC), a 
macroscopic quantum mechanical wave packet with non- 
linear behavior due to interatomic interactions. Such a 
macroscopic coherent matter wave packet of BEC opens 
a new window to study the nonlinear quantum dynam- 
ics governed by the nonlinear Schrodinger equation or 
the mean field Gross-Pitaevskii (GP) equation [lfj. No 
doubt, tunneling of a BEC matter wave packet would 
exhibit different behaviors compared with the single par- 
ticle tunneling in the linear quantum mechanics. In 
fact, there already exist many studies on BEC tunnel- 
ing through different kinds of potentials in the literature 
ll-[l3|. However, to our knowledge, the effect of non- 
linear interparticle interaction on the Hartman effect has 
not been explored. 

In this brief report, we theoretically investigate how 
the interatomic interaction affects the tunneling time 
in the case of a coherent BEC wave packet traversing 
through a potential barrier created by laser beam. We 
consider a BEC wave packet confined in a quasi-one- 
dimensional atomic waveguide to traverse a potential 
barrier created by a far blue-detuned laser beam, which is 
shown in Fig. [TJ Suppose that the quasi-one-dimensional 
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FIG. 1: (Color online) Schematic diagram of a BEC wave packet 
traversing a blue detuning laser beam in an atomic waveguide. 



atomic waveguide is transversely confined by a harmonic 
trapping potential, and the dynamics of atomic BEC can 
be described by the 3D nonlinear Schrodinger equation 
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where 'F is the normalized macroscopic wave function 
of Bose Einstein condensates, rj_ = (y,z), m the 
atomic mass, and co± the trapping frequency in the ra- 
dial(transverse) direction, g = Airh 2 a s /m describes the 
interatomic interaction with a s being the atomic s-wave 
scattering length. TV is the total number of atoms in the 
condensate, and V (x) is the potential barrier created by 
the blue-detuned laser beam. 

When transverse confinement is very strong, the trans- 
verse motion of BEC atoms may be considered to remain 
in the ground state. As a result, we can approximate the 
total wave function of the BEC as 
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where ip± (r) is the transverse ground state wave function 
of the BEC, which can approximately be replaced by a 
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Gaussian function in the weak nonlinear limit 

with a ± = \Jh/mu)±_ being the ground state length 
of harmonic trapping potential. The longitudinal wave 
function ip (x, t) is governed by the effective ID nonlin- 
ear Schodinger equation: 

i^ = -\^ + V{x)iP + g\iP\ 2 ip. (4) 

In Eq. (j4]), we express x in units of a ± , t in units of wj 1 , 
the dimensionless nonlinear interaction g = 2a s N/a ± , 
and the dimensionless wave function ip = ip/^a ± . V (x) 
is the potential barrier experienced by the BEC wave 
packet and written in the form 

V (x) = Vof (x) , (5) 

where Vq is the peak value of the potential normalized by 
Huj± and / (x) is the barrier profile, which is controlled 
by the laser beam. 

Now we have obtained a general dimensionless equa- 
tion dU) to describe the BEC wave packet transmit- 
ting through a barrier. As pointed out in Ref. [lij . 
Eq. ([!]) is invariant under the scaling transformation 
{x, t, ip, Vq} —> {xi], trj 2 ,ip/rj, Vb/?7 2 } where r\ is a dimen- 
sionless constant. To cover different experimental situ- 
ations, it is convenient to introduce such scaling trans- 
formation in numerical simulation. In this paper, we set 
lu±_ R3 2n x 100 Hz and rj = 10, then for 87 Rb atoms, 
a x s» 1 /im. In principle, Eq. (|4|) can describe a BEC 
wave packet transmitting through an arbitrary barrier. 
To simplify the problem, without loss of physical feature, 
in this paper we just consider a rectangular barrier case. 
Such a rectangular barrier can be approximately created 
by a su per- Gaussian laser beam with a large enough or- 
der |15l - [l7| . Therefore the barrier profile / (x) has the 
form 

/(aO = ln ~^l X< ^\ (6) 
J w \ 0, x < — I or x > j. w 

In the following, we simulate finite wave packets 
traversing through the rectangular barrier via the split 
operator method (2lj . Assume the normalized initial 
wave packet is Gaussian 

^ (x ' 0) = exp (" ^aS^ + lko {x ~ xo) ) ' 

(7) 

where xq, Ax and ko are the initial center position, the 
initial half width and the initial center momentum of 
the wave packet, respectively. From Eq. ^ we can 
find that the interatomic interaction will induce the self- 
phase modulation(SPM) [3 [l9[ of matter wave packet. 
The SPM happens in both free region and potential re- 
gion and causes confusion of the nonlinear effect on the 



tunneling process inside the barrier with the SPM pro- 
cess outside the barrier region. To avoid this confusion, 
we shall eliminate the nonlinear interaction outside of 
the potential region, i.e., make the s-wave scattering 
length a s vanishing outside of the potential region. In 
principle, this can be realized by employing the Fesh- 
bach resonance technique to tune the atomic scattering 
length via a spatially varying magnetic field 20]. For 
spatially varying magnetic field, the scattering length 
a s (x) — a s o (1 — AB/(B(x) — Bo)), where a sQ is the 
background scattering length, AB the resonance width, 
Bq the magnetic field of resonance, B(x) the magnetic 
field. By tuning the magnetic field, one can control the 
atomic nonlinearity with positive sign (corresponding to 
repulsive interatomic interaction), negative one (attrac- 
tive interaction) and zero (no interaction). In our work 
we set the spatial profile of the scattering length is con- 
sistent with that of the barrier, namely, g (x) — gof(x). 
Considering the finite width of initial atomic BEC wave 
packet used in experiment, we set Aa; = 50 in all the 
following simulations, corresponding to a half-width 500 
/im for a BEC wave packet composed of 87 Rb atoms. 
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FIG. 2: (Color online) The normalized transmitted wave packets 
(solid curve) at time t = 440 with nonlinear interaction: (a) go = 
— 16, (b)go = and (c)go = 200. Other parameters are Vq = 1, 
L = 12 and ko = 1.2. The freely propagating reference wave packet 
is represented by the dashed curve. 

With the above assumptions, we now investigate the 
tunneling of the BEC wave packet through the barrier 
with atomic nonlinearity. In Fig. 2, we present the trans- 
mitted wave packets with different nonlinear interaction 
at time t — 440, at which the transmitted wave pack- 
ets just emerges from the barrier. Meanwhile we also 
show the freely propagating wave packet as reference, for 
which both the barrier and nonlinear interaction are set 
to zero. From the top panel, one observes that the trans- 
mitted wave packet with negative nonlinearity slightly 
lags behind the reference one. The middle panel is the 
linear case, where the nonlinear interaction is neglected. 
In contrast to the negative nonlinear case, the transmit- 
ted wave packet is ahead of the reference one due to the 
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finite- width effect of incident wave packet, as has been 
pointed out by previous works |22 - [24| . In the bottom 
panel, we find the transmitted wave packet with posi- 
tive nonlinearity is far ahead of the reference one, which 
is quite different from other two cases. The results in 
Fig. 2 show that the nonlinear interaction indeed, as ex- 
cepted, affects the tunneling process of a BEC matter 
wave packet through the barrier. Below we explain the 
physics behind the numerical results in detail. 




FIG. 3: The corresponding momentum spectra of transmitted 
wave packets in Fig. [2] Solid curve stands for the case without 
nonlinearity, and dashed and dash-dotted curves for the cases with 
positive and negative nonlinearity, respectively. 

We consider the transmitted spectra in the momentum 
space. Due to the SPM-induced spectral broadening as 
well as the filtering effect of the barrier 0, , the central 
momenta in the transmitted spectra exhibit great differ- 
ence in the linear, positive and negative nonlinear cases, 
as shown in Fig. [3l The filtering of the barrier takes effect 
on both the linear and nonlinear cases, while the spectral 
nonlinear broadening only occurs in the nonlinear cases. 
The combination of filtering effect and nonlinear broad- 
ening leads to a large transmitted central momentum in 
the positive nonlinear case. Therefore, with the same 
barrier parameters and incident wave number, the BEC 
wave packet with repulsive interatomic interaction tends 
to spend less time in the barrier. However, we note that 
the transmitted central momenta in the negative non- 
linear and linear cases are less than that in the positive 
nonlinear case. This can also be well understood from 
the view of point below. In fact, when the interatomic 
interaction exists, the atomic nonlinearity modifies the 
potential barrier and the atoms " see" an effective poten- 
tial 



V eff = V f {x) + g{x)\i>{x,t)Y 



(8) 



For a repulsive interaction, the nonlinear term has a pos- 
itive sign. Consequently, the height of the linear barrier 
is raised, and vice versa for an attractive interaction. In 
terms of the Hartman's calculation [Ej], the wave packet 
traverses faster through a higher barrier than a lower one 
with the same and enough wide barrier width. Therefore 
the transmitted wave packet with negative or positive 



nonlinearity spend more or less time than the linear one 
in the barrier. 

Now we turn to the quantitative calculation for the 
tunneling time of BEC wave packet through the laser- 
induced barrier. For a finite wave packet transmitting 
through a barrier, we can not use the usual stationary 
phase method to find the transmission time of the trans- 
mitted wave packet, because the condition for the sta- 
tionary phase method is not satisfied in this situation. 
However, the widely used time-of-flight method is appli- 
cable to measure the transmission time of transmitted 
wave packet (26J. The method is described as follows. 
Assume that the incident wave packet is placed at the 
position x (0) = —xq{xq > 0) at the time t = 0, some- 
where to the left of the barrier, and let it move to the right 
with initial momentum fco. After transmitting through 
the barrier the position of transmitted wave packet is 
located in x (ty) at the time t = tx- Due to the ef- 
fect of barrier and nonlinear interaction, the transmit- 
ted wave packet is usually deformed. In this case, the 
appropriate way to describe the center of the transmit- 
ted wave packet is to define the expected position as 



x = I X>0 X \^[ X ^)\ dx/ J x>Q \tfj(x,t)\ z dx [25 
one gets the tunneling time 



So 



M = t T - X °- L l 2 - X ^l- L '\ (9) 



where fco is the transmitted central momentum, defined 
as 



fco = / k\i; T (k)\ 2 dk, 



(10) 



with ipx (fc) the momentum distribution of transmitted 
wave packet. 
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FIG. 4: Tunneling time as a function of the nonlinear interaction 
strength. Other parameters arc fco = 0.6, Vfo = 1 and L = 6. 

First, we study the dependence of tunneling time on 
the nonlinear interaction. Fig. [4] plots the behavior of 
tunneling time via nonlinear interaction strength. Ob- 
viously the tunneling time descends with the nonlin- 
ear interaction strength increasing. This means that 
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the stronger the repulsive interatomic interaction is, the 
faster the BEC wave packet traverses through the barrier. 
It is interesting to note that the tunneling time is negative 
for large positive nonlinearity. The negative tunneling 
time implies that the transmitted wave packet exits from 
the barrier just before the incident wave packet arrives at 
the barrier. Such a phenomenon has also been found and 
discussed in the literatures for different situations [28l |29| . 
Here the negative tunneling time is due to the repul- 
sive interatomic interaction, and completely a nonlinear 
quantum mechanical phenomenon for an atomic Bosc- 
Einstcin condensate. 
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FIG. 5: Tunneling time as a function of the barrier width with 
nonlinear strength go = (Solid line), go = 5 (dashed line) and 
go = —2 (dash-dotted line). Other parameters are fco = 0-6 and 
V = 1. 



Now, we examine the dependence of tunneling time 
on the barrier width as shown in Fig. [5] We find that 
the tunneling time in linear and positive nonlinear cases 
decreases, while increases in negative nonlinear case, with 
the barrier width increasing. These results imply that the 
Hartman effect is violated in both nonlinear and linear 
cases. For the linear case, such a violation is due to the 
finite width of the incident wave packet. This has been 
touched in many theoretical and experimental works [2^- 

In summary, we have numerically studied the effect of 
nonlinear interatomic interaction on the quantum tunnel- 
ing time of a BEC matter wave packet through a laser- 
induced barrier. Analysis shows that both the sign and 
the strength of nonlinear interaction significantly affect 
the tunneling time. As a result, the so-called Hartman 
effect in the linear quantum mechanics could be violated 
in the nonlinear quantum mechanics with a macroscopic 
matter wave packet of Bose-Einstein condensate. 
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